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ABSTRACT 



Chen (1996) has recently shown how the response matrix of a two-phase composite 
can be written as certain linear combinations of products of the component matrices. 
We elaborate on Chen's expansions by deriving them in a different way, which a. 
shows them in a different light, and b. permits us to generalize them further. As an 
application of our results we find exact microstructure-independent relations between 
the moduli of the two components and those of any composite. The body of these 



relations is equivalent to the compatibility relations of Milgrom and Shtrikman (1989a 
ON ' but they are cast in a rather different form, which has certain advantages. As an 

example, we show how any modulus of an arbitrary two-phase composite can be written 
in closed form as a linear combination of any other n of its moduli, with coefficients 
that depend only on the component moduli, but not on the volume fractions, or the 
ON . microstructure. 
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Our discussion concerns the response of a medium to many coupled fields, such as in the 
thermoelectric phenomenon, in the magneto-electric effect, or in coupled, multispecies diffusion. 
We concentrate of the response properties of a composite-whose multifield linear-response 
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matrix is L*-made of two components with response matrices L a , Lb- Chen (1996) , using a 
matrix-decomposition scheme, has recently shown that one can write e.g. 

L* = aoLb + a\L a + aiL a L h L a + ... + a n -\L a L h 1 L a ...L b L a , (1) 

with L a appearing n — 1 times in the last term, were n is the number of coupled fields involved. 
Chen derives several other similar expansions, and there are similar expansions for (L*) . The 
expansion coefficients are functions of the detailed makings of the composite-viz. the volume 
fractions, and the microstructure. 

We rederive and generalize Chen's expansions in §2, and determine their coefficients in §3. In 



4 we present a new, determinant form for the compatibility relations of Milgrom and Shtrikman 
(1989a). 
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2. Rederivation and generalization of Chen's expansions 



Following Milgrom and Shtrikman (1989a) we take W to be a (non-singular) matrix that 
diagonalize L a and L b simultaniously by a congruence transformation, such that 

WL a W = A, WL b W = I, (2) 

where A = diag(\i, A n ) is diagonal, and / is the unit matrix. Such a matrix always exists as 
L a , L b are symmetric and positive definite. As explained in Milgrom and Shtrikman (1989a), the 
physics of the problem dictates that the response matrix, L* , of any composite material of the two 
components is also digonalized by W: 



WL*W = A* = diag{\\, A*). 



(3) 



Consider now the series ofnxn matrices P", for k an arbitrary integer, such that is 
transformed by W into A fc -the A:th power of A: 



P (k) 



(4) 



It is easy to see from eq.@ that P^ = L b , P^ = L a , and, in general, 

P( fc ) = L b (L^L a ) k . 



(5) 



Thus, for positive k P^ = L a L h l L a ...L b l L a , and for negative k = —\k\ P^ = LbL a 1 Lf ) .-.L a 1 Li ) , 
with L a , or L^ 1 appearing \k\ times. If one interchanges the role of L a and L b (W is than 
changed): P^(L a , L b ) = P^ k+1 \L b , L a ). 

Let N be the number of different AjS, than it is clear that: 

a. Any set of more than N integer powers of the matrix A are linearly dependent (because 
each can be viewed as a vector of length N). 

b. Any set of different (integer) powers, A fcl ,...,A fcjv are linearly independent. This is 
because the determinant of the matrix 



M(A, k) = 



I Ai 1 



A^\ 



\fcjv / 



(6) 



does not vanish unless at least two of the Aj are equal (for the case of consecutive powers this is 
just proportional to the Vandermonde determinant-see below) (see e.g. Aitken 1949 )| ) . 



c. Any n x n diagonal matrix, with the same equalities among its element as in A, can be 
written as a linear combination of any set of N different powers of A. Property c. is satisfied, in 
particular, by A* for any composite, by its inverse, or, for that matter, by any integer power of it. 
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For the sake of clarity we assume from now on that A is non-degenerate i.e. N = n; all we say 
below goes though straightforwardly for N < n. 

Properties a.-c. are carried through by transforming back with W, whereby X k goes back to 
P<*\ /(A) is carried back to f(L a L^)L b = L b f(L^L a ), X* goes to L*, (A*)" 1 goes to L b {L*)- l L b: 
(A*) m (to positive) goes to L* L^ 1 L* ...L^ 1 L* (L* appearing m times), etc. We thus have 

a'. Any set of more than n of the P^ matrices are linearly dependent. 

b'. Any set of n different p( k 's are linearly independent. 

c'. Any of the matrices, such as L*, L b (L*) L b , etc. can be written as a linear combination 
of any n of the P( k >s (any N when N < n). 

So, for example, for every choice of different powers there are coefficients ai such that 

n 

L* = Y j a i P {kl) . (7) 

i=l 

Chen's expansions are special cases of eq.(^) [and an analogous one for L b (L*)~ 1 L b ] with 
consecutive values of kf. k = ko, ko + 1, ko + n — 1, and special values of ko. For example, eq.(|l]) 
is obtained when expanding L* with ko = 0. Expanding L^L*)^ 1 Lj, using ko = —n + 1 gives 

(L*)- 1 = fto^ 1 + hL- 1 + b 2 L- 1 L b L~ 1 + ... + b n ^L- l L b ...L b L-\ (8) 

another of Chen's expansions. 

Obviously, instead of using the P( k h as expansion basis, we could use the response matrices 
of different composites of the same two composites, or different space-space moduli matrices 
of anisotropic composites of the same components. This will give relations between different 
composites [for more details see the discussion in Milgrom and Shtrikman (1989a)]. 



3. The expansion coefficients 



The expansion coefficients aj depend on the exact nature of the composite. It is sometimes 



useful to express these coefficients in terms of Aj and A*, which |Chen (1996)| does using a matrix 
decomposition theorem. In the present formalism the ctj are just the coefficients of the linear 
expansion of the n- vector A*, made of the diagonal of A*, in the n n- vectors made of the diagonals 
of \ ki . The vector a = (ai, a n ) is thus the solution of the linear equations 



where M is the matrix defined in eq. 
1949)1 (p. Ill)] 



M{X,k)a = X*, (9) 
[Their ilk are treated under 'alternant matrices' in Aitken| 
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The case of consecutive ks lands itself to further simplification: When 
k = (ko, ko + 1, &o + n — 1) eq.(P) can be written as 



where V is the Vandermonde matrix: 



V(X)a 



/l Ai 
1 A 2 



(*o) 



F(A) 



A? 
A| 



AT 1 \ 

\ n— 1 
A 2 



and 5* = A*/A, fco . Thus 



V 1 A n A^ 
4* ») = £(y-%-A*/(A 



AIT 1 / 



The elements of V 1 can be written in closed form: 

= (-l) i+J '|^n|- 1 |^- 1 |S , „-i(Ai,...,A J -_i,A i+1 , 



■j A„). 



(10) 



(11) 



(12) 



(13) 



Here, \V n \ is the Vandermonde determinant for the n Aj, |V^_ 1 | is the same for the n — 1 As 
excluding Aj, and S n ~i is the elementary symmetric polynomial of degree n — i in its variables 
(here Ai, Aj_i, Aj + i, A n ); it is simply the sum of all possible products of n — i subsets of 
different variables. So, for example, Sq = 1, S± = J2k^j Afc> and 5 n _i = AiA2...Aj_iAj + i...A n . All 
this gives the expressions Chen (1996) obtains for the special values of ko corresponding to his 
expansions. 



4. A new form of the compatibility relations 



The fact underlying many an exact relation for two-phase composites is that only n numbers 
are needed to specify all the the n(n + l)/2 moduli of a response matrix. Thus, the off-diagonal 
elements of eq.(||) can be viewed as n(n — l)/2 linear relations in the elements of L* with 
coefficients involving only the component moduli (through W). A given set of the compatibility 



relations of Milgrom and Shtrikman (1989a) is another, and more useful, set of n(n — l)/2 relations 
between the moduli. In Chen's expansions and all their generalization discussed above [eq.(0)] the 
n structure-dependent parameters are the expansion coefficients. While all these sets of constraints 
are equivalent, each has its advantages. We now proceed to derive yet other manifestations of the 
same constraints, which have their own advantages. 

Let A A n+l be any n + 1 matrices that become diagonal when transformed with W . They 
can be any of the P( k 's, the response matrices of any isotropic composite of the two components, 
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the space-space response matrix of an anisotropic composite, or appropriate products of such 
matrices of the form Mi M^ 1 M%...M^_ 1 M m , where the Ms are diagonalized by W. The 
arguments of §1 show that these n + 1 matrices are linearly dependent. Thus, if we choose any set of 
re + 1 matrix slots: risi, r n +is n+ %, the re + 1 vectors a, = [-AJ. iai , A%- n+lSn+1 ]i (1 < i < n. + 1), 
are linearly dependent: 

A 1 



A 1 



A n+1 



An+1 

Vn+lSn+l 



0. 



Take, for example, for the above n + 1 matrices any n different P^s with k = k±, 
response matrix, L*, of some isotropic composite. Then eq.(14) tells us that 



(14) 



.,k n , and the 



T* 



P 



rxsi 



T* 

r 2 S2 



(fci) p(fci) 
it 



T* 

r n+1 s n+ i 



P 



(*i) 



0, 



(15) 



for any choice of different matrix slots rjSj. 

Equation (|l~5|) constitute a linear relation, with structure-independent coefficients, between 
any chosen set of n + 1 composite moduli. Any set of independent n(n — l)/2 such determinant 
relations is equivalent to one matrix compatibility relation of |Milgrom and Shtrikman (1989a) of 
the form 

VL- a x L h = L h L- a x L\ (16) 

which embodies the same number of independent relations between the elements of L a , and L*. 
The advantage of the latter is that it involves only elements of the products L~ 1 L; ) , and L^L" 1 , 
while the former involves higher-order products P^ . The advantage of the former is that each of 
its relations involve only n + 1 of the elements of L*, and they can be chosen at will, while in the 
latter each relation involves generically 2n — 1 elements of L*, which cannot be chosen at will. The 
new relations eq.(|l5|) thus permit us to write any modulus of L* as a microstructure-independent 
linear combination of any other n of its moduli. (Clearly, if we replace L* in eq.(|l5|) by another of 
the P( k h we get an identity, as the derivation of such a relation involves no physical arguments.) 

One interesting example is the choice with n of the slots being the diagonal ones, and the 
(n + l)th some off-diagonal element rs. This gives a determinant relation of the form 



r * /* 
p(fel) _(fcl) 



rs Pi 



1* 

D (*i) 
Pn 



„(*») 
Pn 



0, 



(17) 
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where I* are the diagonal coefficients of L* , and pi is the jth diagonal element of P^ ki \ This 
form permits us to express all the off-diagonal moduli of any composite in terms of the diagonal 
moduli (and those of the components): 

L* rs = J2<x r i S lh where <^' = ^-. (18) 

(k) (rs) 

Here, \P\ is the determinant of the diagonal elements p, , and 17} are the appropriate cofactors, 
and all are calculated from the component moduli. An off-diagonal modulus of an arbitrary 
composite is then given by a universal linear function of the diagonal moduli of this composite. 



For the special case when only two coupled fields are involved Milgrom and Shtrikman (1989a) 
were, in fact, able to write the single compatibility relation in determinant form (in this case 
2n — 1 = n + 1). Our present result is a generalization to the n > 2 case. 

Equation© can be useful, for instance, in the following way: If all cross moduli are 
small and can be neglected to lowest order, one can find simultaneous bounds on the diagonal 



elements I* pa A*: A region B in the space of ■■■,1%)- Equation(18), with the left-hand side a 
parameter, defines a family of hyperplanes in this space. Bounds on individual L* s can then be 
found by identifying the two extreme hyperplanes of the family still touching the region B. Note 
that these would constitute (approximate) bounds on individual cross moduli. In contrast, the 
bounds discussed by Milgrom and Shtrikman (1989c) , and by |Chen (1996) , are exact but involve 



combinations of many moduli. 

All our results are applicable also to polycrystals made of a uniaxial single crystal as is evident 



from the discussion of Milgrom and Shtrikman (1989b) 



I thank Tungyang Chen for permission to use and quote his results before publication. 
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